Nonlinear dynamics and nonlinear dynamical systems 

Controlled Hopf Bifurcation of a storage-ring free-electron laser 
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Local bifurcation control is a topic of fundamental importance in the field of nonlinear dynamical 
systems. We discuss an original example within the context of storage-ring free-electron laser physics 
by presenting a new model that enables analytical insight into the system dynamics. The transition 
between the stable and the unstable regimes, depending on the temporal overlapping between the 
light stored in the optical cavity and the electrons circulating into the ring, is found to be a Hopf 
bifurcation. A feedback procedure is implemented and shown to provide an effective stabilization 
of the unstable steady state. 

PACS numbers: 05.45.-a 
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Transition from stability to chaos is a common char- 
acteristic of many physical and biological systems 0, Q ■ 
Within this context, local bifurcation control is a topic 
of paramount importance especially for those systems 
in which stability is a crucial issue. This is, for ex- 
ample, the case for conventional and non conventional 
light sources, such as Storage-Ring Free-Electron Lasers 
(SRFELs), commonly employed in various scientific ap- 
plications In a SRFEL the physical mechanism 
responsible for light emission and amplification is the in- 
teraction between a relativistic electron beam and the 
magnetostatic periodic field of an undulator. Due to the 
effect of the magnetic field the electrons emit synchrotron 
radiation, known as spontaneous emission. The light pro- 
duced by the electron beam is stored in an optical cavity 
and amplified during successive turns of the particles in 
the ring. A given temporal detuning, i.e. a difference be- 
tween the electron beam revolution period and the round 
trip of the photons inside the optical cavity, leads to a cu- 
mulative delay between the electrons and the laser pulses: 
the laser intensity may then appear as "continuous wave 
(cw)" (for a weak or strong detuning) or show a stable 

galsed behavior (for an intermediate detuning amount) 
6]. The achievement of a large and stable "cw" zone 
is a crucial issue, of fundamental importance for experi- 
mental applications. In this Letter, we characterize the 
transition between stable the unstable regimes as a Hopf 
bifurcation. This result allows one to establish a formal 
bridge with the field of conventional lasers and to adopt 
the universal techniques of control theory to enlarge the 
stable signal region. We develop this idea by introducing 
a new model which reveals to be particularly suitable for 
analytic investigations. 

The longitudinal dynamics of a SRFEL can be described 
by a system of rate equations accounting for the coupled 
evolution of the electromagnetic field and of the longitu- 
dinal parameters of the electron bunch 0, . The tem- 
poral profile of the laser intensity, is updated at each 



pass, n, inside the optical cavity according to: 

y„+i(r) = i?2y„(r - e) [1 + .g„(r)] + ^.(t), (1) 

where r is the temporal position of the electron bunch 
distribution with respect to the centroid; R is the mir- 
ror reflectivity; the detuning parameter e is the difference 
between the electrons revolution period (divided by the 
number of bunches) and the period of the photons inside 
the cavity; is{T) accounts for the profile of the sponta- 
neous emission of the optical klystron 9] . Assuming that 
the saturation is achieved when the peak gain is equal to 
the cavity losses, P, the FEL gain 5„(t) is given by 0,13: 
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where gi and (Tq are the initial (laser-off) peak gain and 
beam energy spread, (t„ and the energy spread 

and the bunch length after the nth light-electron beam 
interaction, and 7 = — (Jq. Note that equation (0) 
refers to the case of SRFELs implemented on an optical 
klystron. The evolution of the laser-induced electron- 
beam energy spread is governed by the following equa- 
tion: 



2AT, 



(3) 



Here cTg is the equilibrium value (i.e. that reached at the 
laser saturation) of the energy spread at the perfect tun- 
ing and AT is the bouncing period of the laser inside the 
optical cavity; /„ = J_^yn{T)dT is the laser intensity 
normalized to its equilibrium value (i.e. the saturation 
value for e = 0) and stands for the characteristic time 
of the damped oscillation of electrons in their longitudi- 
nal phase-space. Equations JQ), Q and are shown 
to reproduce quantitatively the experimental results 0. 
In particular, the laser intensity displays a stable "cw" 
behavior for small amount of detuning, while a pulsed 
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regime is found for e larger than a certain critical thresh- 
old, Ec- This model represents the starting point of om' 
analysis. 

Equation ^ characterizes the evolution of the statis- 
tical parameters of the laser distribution: by assuming 
a specific form for the profile, it is in principle possible 
to make explicit the evolution of each quantity. For this 
purpose, we put forward the assumption of a Gaussian 
laser profile and compute the first three moments. The 
details of the calculations are given elsewhere jl^- In 
addition, it is shown that for e spanning the central "cw" 
zone, the quantities [pi^nl <^T,n)^ and [(t„ -|- e)/(T7-_„]^ are 
small. Hence, a Taylor series expansion is performed and 
second order terms neglected. Finally, by approximating 
finite differences with differentials, the following contin- 
uous system is found: 



energy spread, intensity, centroid position and rms value 
of the laser distribution, to the light-electron beam de- 
tuning. Through a stability analysis it is also possible to 
determine the threshold value t^. To our knowledge this 
study represents the first attempt to fully characterize 
the detuned SRFEL dynamics. 
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FIG. 1: Phase-space portraits. Left column: I is plotted ver- 
sus I. Right column: a versus a. The top panels refer to 
e — 0.1/s < Ec, the bottom ones to e — 1.3/s > ec. Simu- 
lations have been performed using the case of the Super-ACO 
PEL as reference. The values of the relevant parameters are: 
AT = 120 ns, Ts = 8.5 ms, ao = S-IO"", ae/ao = 1.5, = 14 
kHz, g, = 2%, P = 0.8%, h = 1.4 ■ 10"^ 

The fixed points (/, ct, r, 07) are found by imposing 
f = f = f = IF = in 0, and solving the cor- 
responding system. Assume hereon e > 0, being the 
scenario for e < completely equivalent. After some 
algebraic calculations, the following relations are found: 



Here VI represents the oscillation frequency of the elec- 
trons in their longitudinal phase-space and a, the mo- 
mentum compaction factor, is a characteristic parameter 
of the storage ring. Note the redefinition of a which is 
from hereon normalized to ctq. Although in approximate 
form, system Q still captures the main features of the 
longitudinal SRFEL dynamics. In particular, the transi- 
tion from the "cw" regime to the unstable (pulsed) steady 
state occurs for a temporal detuning which is close to 
the one found in the framework of the exact formulation, 
hence to the experimental value. However, system (0J 
fails in reproducing the correct behavior when the tran- 
sition to the lateral "cw" zone is approached. In Figure 
n phase-space portraits for both the laser intensity and 
the beam energy spread are plotted for different values 
of e. Limit cycles are observed when e > Cc. For smaller 
values of e, the variables converge asymptotically to a 
stable fixed point. The latter can be analytically char- 
acterized, thus allowing one to relate the electron-beam 
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These relations link the equilibrium values of I,T,ai to 
a. The quantity ct is found from the following implicit 
equation: 



Laser intensity 



Electrons' energy spread 
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where o7 and r are respectively given by © and ((SJ. 
For any given value of the detuning e, equation Ijllll 
can be solved numerically, by using a standard bisection 
method. The estimates of a are then inserted in equa- 
tions (71, to compute the corresponding values 
of /, r, 57. Results of the calculations (solid line) and 
direct numerical simulations using the system Q (sym- 
bols) are compared in Figure 12 displaying remarkably 
good agreement. It is worth stressing that, by means of 
a perturbative analysis, a closed analytical expression for 

as a function of e is also found. The details of the quite 
cumbersome calculations are given elsewhere (Id] . 

As a validation of the preceding analysis, we consider 
the case of perfect tuning, i.e. e = 0, and compare our 
estimate for the laser induced energy spread 57 to the 
value (57) sm, derived in the context of the widely used 
super-modes approach Both theoretical predictions 
are then compared to experiments performed on the Su- 
perACO and Elettra FELs. Results are given in TableQ] 
the improvement of the calculation based on equation (jSJ 
is clearly shown. The results of Tablenjindicate that both 
07 and {Wi)sm are smaller than the experimental values. 
This is probably due to the fact that the models neglect 
the effect of the microwave instability ^2 resulting from 
the electron beam interaction with the ring environment 
(e.g. the metallic wall of the vacuum chamber). In the 
case of Elettra the situation is complicated by the pres- 
ence of a "kick-like" instability (having a characteristic 
frequency of 50 Hz) which periodically switches off the 
laser preventing the attainment of a stable "cw" regime 

la. 





Super AGO 


Elettra 


m (ps) 


5 


2 


(o7)sm (ps) 


3 


1 


Experimental values (ps) 


10 ±2 


5 ± 2 



TABLE I: Theoretical widths of the laser pulse compared to 
experimental values for the case of the Super-ACO and Elettra 
FELs. The experimental setting for the case of Super AGO 
(operated at a beam energy of 800 MeV and at a laser wave- 
length of 350 nm) is that specified in the caption of Figure 
The analogous parameters for ELETTRA (operated at a 
beam energy of 900 MeV and at a laser wavelength of 250 nm) 
are the following: AT — 216 ns, Ts — 87 ms, ao = 1 • 10""^, 
(Je/cro = 1.5, Q. = m kHz, Qi = 15%, F = 7%, Is = 4.3- 10"^ 



The stability of the fixed point [/(e), o'(e),T(e), 07(e)] 
can be determined by studying the eigenvalues of the Ja- 
cobian matrix associated with the system Q. The real 
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FIG. 2: The fixed points are plotted as function of the de- 
tuning parameter e. Top left panel: Normalized laser in- 
tensity. Top right panel: Normalized electron-beam energy 
spread. Bottom left panel: Laser centroid. Bottom right 
panel: rms value of the laser distribution. Symbols refer to 
the simulations, while the solid line stands for the analytic 
calculation. The list of parameters is enclosed in the caption 
of Figure^ 



part of the eigenvalues as a function of e is shown in Fig- 
ure 121 The system is by definition stable when all the 
real parts are negative. The transition to an unstable 
regime occurs when at least one of them becomes posi- 
tive. In general, the loss of stability takes place according 
to different modalities. Consider the case of a Jacobian 
matrix with a pair of complex conjugate eigenvalues and 
assume the real parts of all the eigenvalues to be negative. 
A Hopf bifurcation occurs when the real part of the two 
complex eigenvalues becomes positive, provided the other 
keep their signs unchanged 14]. This situation is clearly 
displayed in Figure O thus allowing to conclude that the 
transition between the "cw" and the pulsed regime in a 
SRFEL is a Hopf bifurcation. The critical detuning, Ec, 
can be calculated (open circle in Figure Q) and displays 
good agreement with both the simulated data and the 
experimental value. A closed relation for Ec is also found 
[Toi, by making use of the analytic expressions for the 
fixed points. 

Having characterized the transition from the stable to 
the unstable steady state in term of Hopf bifurcation 
opens up interesting perspectives to stabilize the signal 
and dramatically improve the system performance. In or- 
der to maintain the laser-electron beam synchronism and 
avoid the migration towards one of the unstable pulsed 
zones of the detuning curve, existing seco nd-g eneration 
SRFELs, such like Super-ACO and UVSOR have 
implemented dedicated control systems. The idea is to 
re-adjust periodically the radio-frequency, thus dynami- 
cally confining the laser in the central "cw" zone. Even 
though generally suitable for second-generation SRFELs, 
these systems are inappropriate for more recent devices, 
such as ELETTRA and DUKE. The latters are indeed 
characterized by a much narrower region of stable signal 
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FIG. 3: Real part of the eigenvalues of the Jacobian ma- 
trix associated to the system ^ as a function of the detuning 
parameter e. The solid line refers to the complex conjugate 
eigenvalues. The circle represents the transition from the sta- 
ble to the pulsed regime, i. e. the Hopf bifurcation. 

only occasionally experimentally observed , making a 
priori impossible to pursuit the former strategy. On the 
contrary, the approach here discussed exploits an univer- 
sal property of SRFELs, thus allowing to overcome the 
limitations of other schemes. The procedure consists in 
introducing a specific self-controlled (closed loop) feed- 
back to suppress locally the Hopf bifurcation and enlarge 
the zone of stable signal. This is achieved by replacing 
the constant detuning with the time-dependent quantity 

m- 

e{t) = eo + PATi , (12) 

which is added to system Q. Here eo is assumed to be 
larger that e^. when the control is switched off, i.e. (3 = 0, 
the laser is unstable and displays periodic oscillations. 
For (3 larger than a certain threshold, /3c, the oscillations 
are damped and the laser behaves as if it were in the "cw" 
region. Note that, as soon as saturation is reached, 1 = 
and, thus, the stable regime is maintained asymptotically 
for e = eo > ec, i.e. well inside the former unstable zone. 
The results of the simulations are represented in Figure 
il 

This new theoretical insight sets the ground for exper- 



imental tests [lO|. In this respect, a significant and re- 
producible extension of the stable "cw" region using this 
technique has been recently achieved at Super AGO 0|. 
This result fully confirms our theoretical predictions. 

In conclusion, in this Letter we propose a new approxi- 
mate model of a SRFEL. This formulation enables a deep 
analytical insight into the system dynamics, allowing one 
to derive the explicit dependence of the main laser pa- 
rameters on the temporal detuning. Results are fully con- 
firmed by numerical simulations and show satisfactory 
agreement with available experimental data. Further, 
the transition between the stable and unstable regimes 
is found to be a Hopf bifurcation, and the critical de- 
tuning Ec is calculated explicitly. Finally, we introduced 
in the model a derivative feedback that is shown to sta- 
bilize the laser intensity well beyond the threshold Ec- 
Successful experiments carried out at Super AGO con- 




FIG. 4: Behavior of the FEL (normalized) intensity in 
absence (upper panel) and in presence (lower panel) of the 
derivative control system. The simulations refer to the case 
of Super AGO (see caption of Figure^ for the list of the pa- 
rameters). Here eo — 1.3 fs > ec. The stabilization has been 
achieved using /9 = 6 ■ 10""^. Here, jSc — 5 ■ 10~*. 



firmed our predictions. Preliminary experiments carried 
out at ELETTRA have also given encouraging results. 
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